Thermal expansion and Griineisen parameter in quantum Griffiths phases 
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We consider the behavior of the Griineisen parameter, the ratio between thermal expansion and 
specific heat, at pressure-tuned infinite- randomness quantum-critical points and in the associated 
quantum Griffiths phases. We find that the Griineisen parameter diverges as ln(Tb /T) with vanishing 
temperature T in the quantum Griffiths phases. At the infinite-randomness critical point itself, the 
Griineisen parameter behaves as [ln(To/7 1 )] 1+1 where v and are the correlation length and 
tunneling exponents. Analogous results hold for the magnetocaloric effect at magnetic-field tuned 
transitions. We contrast clean and dirty systems, we discuss subtle differences between fsing and 
Heisenberg symmetries, and we relate our findings to recent experiments on CePdi-^Rh^. 

PACS numbers: 71.27.+a, 75.10.Nr, 75.40.Cx, 71.10.Hf 
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Introduction — Quantum phase transitions (QPTs) oc- 
cur at zero temperature when a parameter such as pres- 
sure or magnetic field is varied. At these transitions, the 
quantum fluctuations associated with the competition 
between different quantum ground states lead to uncon- 
ventional thermodynamic and transport properties. 1,2 ' 3 ' 4 
In metallic systems they can induce, e.g., non- Fermi liq- 
uid behavior and exotic superconductivity. The charac- 
terization of QPTs is a topic of great current interest 
with many fundamental questions remaining unresolved. 

Over the last few years, the Griineisen parameter T, 
the ratio between thermal expansion coefficient and spe- 
cific heat, has become a valuable tool for analyzing 
pressure-tuned QPTs. For transitions tuned by magnetic 
field, the same role is played by the magnetocaloric ef- 
fect. Zhu et al. 5 ' 6 showed that the thermal expansion 
coefficient is more singular than the specific heat at a 
generic clean quantum critical point (QCP). They thus 
predicted that the Griineisen parameter diverges when 
approaching criticality. Specifically, if hyperscaling holds 
(below the upper critical dimension), F ~ T _1 /( Ziy ) if the 
temperature T is lowered at the critical pressure p c and 
r ~ l/(p— Pc) if the pressure p approaches p c at zero tem- 
perature (z denotes the dynamical exponent). Above the 
upper critical dimension, T still diverges, but the func- 
tional form is modified by dangerously irrelevant vari- 
ables. Diverging Griineisen parameters have since been 
observed at several 7,8 ' 9 ' 10 magnetic QCPs. 

Since many materials feature considerable amounts of 
quenched randomness, the study of QPTs in random 
systems has received much attention recently. The in- 
terplay between quantum fluctuations and static ran- 
dom fluctuations results in more dramatic disorder effects 
at QPTs than at classical transitions, including quan- 
tum Griffiths singularities, 11 ' 12,13,14,15 activated dynami- 
cal scaling, 16,17,18 and smeared transitions. 19,20 A review 
of some of these phenomena can be found in Rcf. 21. In 
view of the insight about the character of a QPT that 
can be gained from the Griineisen parameter and the 
magnetocaloric effect, it is desirable to determine their 
behavior within these unconventional scenarios. This is 



particularly timely because exotic scaling behavior com- 
patible with many predictions of the quantum Griffiths 
scenario has recently been observed 22,23 at the ferromag- 
netic QPT in CePdi-zRh^. 

In this paper, we therefore calculate the thermal ex- 
pansion coefficient and the Griineisen parameter (for 
pressure-tuned transitions) as well as the magnetocaloric 
effect (for magnetic-field tuned transitions) at infinite- 
randomness QCPs and in the associated quantum Grif- 
fiths phases. We use two methods, a heuristic rare region 
theory and a scaling analysis of the QCP itself. 

We define the Griineisen parameter 5,24 T as the ratio 
between the thermal volume expansion coefficient 



13 = V- 1 {dV/dT) p = -V- 1 (dS/dp) q 
and the molar specific heat 

c p = TN- 1 (dS/dT) . 



(1) 



(2) 

Here, V is the volume, N is the particle number, and S 
denotes the entropy. Thus, 

B _ (dS/d P ) T ( ] 

c p V m T{dS/dT) p 1 > 

with V m = V/N the molar volume. For a pressure-tuned 
transition, (dS/dp)T = p^ 1 (95 '/ 'dr)r explores the de- 
pendence of the entropy on the dimensionless distance 
from criticality, r = (p — p c )/p c . For a transition tuned 
by magnetic field H with r = (H — H c )/H c , the same 
dependence is encoded in (dS/dH)? = (dM / dT)fj with 
M the total magnetization. We thus define the magnetic 
analog of the Griineisen parameter, 

(dM/dT) H _ {dS/dH) T _ 1 /<9T' 



T H = - 



T(dS/dT) H T \dH 



(4) 



which can be determined from the magnetocaloric effect. 

Rare region theory — For definiteness, we consider a d- 
dimensional quantum Landau-Ginzburg- Wilson (LGW) 
free energy functional for an ri-componcnt order param- 
eter field 4>. The action of the clean system is given by 25 



S 



dx dy 4>(x) K(x,y) (f>(y) + u I dx </> 4 (x) . (5) 



2 



quantum critical ' 

\ n ,t m,iti/vii' quantum disordered 
\+[ln(Vr)] / ^ 1 

const. -ln(r /I)\ /+ln(7' /7') const. 



ordered phase | 



SO 



WO 




QCP 



WD 



SD 



FIG. 1: (Color online:) Schematic phase diagram close to 
an infinite-randomness QCP. SO and SD denote the strongly 
ordered and disordered bulk phases while WO and WD are 
the weakly ordered and disordered quantum Griffiths phases. 
The terms in the figure show the temperature-dependence of F 
or Th (for pressure and field-tuned transitions, respectively). 



Here, x = (x, t) comprises position x and imaginary 

time r, and J dx = J dx J^ T dr. The Fourier trans- 
form of the bare inverse propagator (two-point vertex) 
K(x,y) reads K(q,u n ) = (r + q 2 + 7|w„| 2 / 20 ) with r 
the bare distance from the (clean) QCP. To introduce 
quenched randomness, we dilute the system with non- 
magnetic impurities of spatial density b, i.e., we add a 
potential, 5r(x) — J2i ~~ X WL to r o- Here, x(i) are 
the random positions of the impurities, and V(x) is a 
positive short-ranged impurity potential. 

In our disordered LGW theory, quantum Griffiths 
phases occur in Ising systems (n = 1) with dissipation- 
less dynamics (zq = 1) or in continuous-symmetry sys- 
tems (n > 1) with Ohmic dissipation (z — 2). 21 We first 
focus on the Ising case, minor differences for n > 1 will 
be discussed later. 

We start our analysis with the weakly disordered quan- 
tum Griffiths phase (WD in Fig. 1). Despite the dilution, 
there are large spatial regions devoid of impurities. They 
can be locally in the magnetic phase even though the 
bulk system is still nonmagnetic. The probability w of 
finding such a rare region or cluster of linear size Lrr is 
exponentially small in its volume, w ~ cxp(— bL RR ). Be- 
cause the cluster is locally ordered it acts as a two-level 
system with an energy gap e ~ cxp(— aL RR ). Combin- 
ing the two exponential laws, we obtain the well-known 
power-law density of states (see, e.g., Ref. 21), 



Pi*) 



c A(r)-l 



(6) 



We have parametrized the nonuniversal power law in 
terms of the Griffiths exponent A = b/a. It vanishes at 
the QCP and increases with increasing distance r from 
criticality. To determine the rare-region contribution to 
the entropy at temperature T, we note that each rare re- 
gion with e < T contributes an entropy of In 2 while those 
with e > T do not contribute significantly. Allowing for 



an ?*-dcpendcnt prcfactor, we thus find 

S(r,T) = Ng(r) (T/T ) x ^ 



(7) 



where To is a microscopic temperature scale. 

Thermal expansion coefficient and specific heat can 
now be calculated easily by taking the appropriate deriva- 
tives of the entropy. From (1) and (2), we find the leading 
low-temperature behavior to be 



1 



g(r)X'(r)(T/To) x ^HT /T) , (8) 

(9) 



V m p, 

g(r)X(r) (T/T ) A W , 



where A'(r) denotes the derivative of A with respect to r. 
In the Griincisen ratio, the temperature dependencies of 
(3 and c p almost completely cancel, resulting in 



r = — 



1 A'(r) 
V m Pc A(r) 



ln(T /T) 



(10) 



The rare region contribution to the Gruneisen parameter 
diverges logarithmically with decreasing temperature in 
the entire WD quantum Griffiths phase. Because A in- 
creases with r, both the thermal expansion coefficient and 
the Gruneisen parameter are positive. This agrees with 
the notion that the low-temperature entropy decreases 
with increasing distance from criticality. 

In the weakly ordered (WO) quantum Griffiths phase, 
the relevant degrees of freedom are strongly coupled clus- 
ters that are sufficiently isolated from the (ordered) bulk 
system so that they can fluctuate independently. This 
requires that the effective coupling J c g of the cluster to 
the bulk is smaller than its energy gap which still reads 
e ~ exp(— aL RR ). To isolate the cluster, it must thus 
be surrounded by a large spatial region that is locally 
in the nonmagnetic phase. Generically, the correlations 
will drop off exponentially with distance in this region. 
The condition J e g < e thus implies that the linear size 
of the isolating region must vary as ln(l/e) ~ Lf lR with 
the cluster size Lrr. We conclude that the probability of 
finding a sufficiently isolated cluster of size L_ra drops off 
as if ~ exp[—b(L RR ) d ], i.e., much faster than in the WD 
phase. The resulting density of states takes the form 



P(e) 



1 



■ cxp 



-A(r)ln d (e /e) 



(11) 



with a nonuniversal A(r) which is the analog of the Grif- 
fiths exponent A(r). Thus, we still obtain a gapless spec- 
trum, but the singularity is weaker than in the WD phase 
in all dimensions d > 1. In particular, the density of 
states vanishes faster than any power law with e — > 0. We 
emphasize that (11) is the generic result, special types of 
randomness can lead to stronger singularities. For in- 
stance, in a percolation scenario (site or bond dilution 
of a lattice), a shell of empty sites or bonds is suffi- 
cient to completely isolate a cluster. In this case, w ~ 
cxphW^ 1 ] and p(e) ~ e" 1 exp[-A» ln^^^eo/e)] giv- 
ing rise to a singularity even stronger than in the WD 
phase. 26 
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The rare-region contribution to the entropy can be es- 
timated as above by simply counting the clusters with 
energy gap e < T. This gives 



S(r, T) = Ng{r) exp -A(r) ln d (T /T) 



(12) 



Using (1) and (2), we obtain the leading low-temperature 
behavior of thermal expansion and specific heat, 







g(r)X'(r) 



■ exp 



V m p, 

c p = d5(r)A(r)exp[-A(r)ln d (T /T) 



A(r)ln d (T /T) ln d (T /T) , (13) 



ln d - 1 (T /T)(14) 



This results in a Griineisen parameter of 







1 



dV m p c A(r) 



HTo/T) 



(15) 



which is (except for the extra factor 1/d) identical to the 
WD result (10). Note that a different power of ln(eo/e) in 
the exponent of (11) (as discussed above for percolation 
disorder) would not change the temperature-dependence 
of r. In the WO quantum Griffiths phase, A decreases 
with increasing r (approaching the QCP). The thermal 
expansion coefficient and the Griineisen parameter are 
therefore negative, again in agreement with the entropy 
decreasing with increasing distance from criticality. Note 
that the rare region contributions to both (3 and c p in the 
WO phase are only weakly singular for d > 1. Therefore, 
they may be sublcading to contributions from other soft 
modes in the system, making (15) hard to observe exper- 
imentally. 

Scaling analysis — We now complement the heuristic 
rare region theory by a scaling analysis of the Griineisen 
parameter at infinite-randomness QCPs. These exotic 
critical points emerge from the strong-disorder rcnormal- 



lzation group 



16.27,28 



and generally occur in conjunction 



with the quantum Griffiths phases discussed above. 21 

According to the strong-disorder renormalization 
group, the density of independent clusters surviving at 
temperature T scales like an inverse volume and thus 
has the scaling form 29 

n(r, T) = MT Q /T)]- d ^ $ (r [ln(T /T)] 1 /(^) > ) . (16) 



where v and ip are the correlation length and tunneling 
exponents. The scaling function <&(?/) is analytic at y > 
(because there is no finite-temperature phase transition 
at r > 0). For small y, we can thus expand Q(y) = $(0) + 
y<£'(0) + . . .. For large positive y (in the WD quantum 
Griffiths phase), $(y) = Ay dv expi-cx^) with A and c 
constants. The scaling function <&(?/) has a singularity 
at some y c < marking the transition to the ordered 
phase. This immediately gives the unusual form of the 
phase boundary, T c (r) = To cxp[— (y c /f) ), sketched in 
Fig. 1. Since the surviving clusters are essentially free, 
each contributes so = In 2 to the entropy. The scaling 
part of the entropy thus reads 

S(r,T) = Ns {HTo/T)}- d/xl '<S> (r MTo/T)] 1 ^) (17) 



We first calculate the thermal expansion coefficient and 
the specific heat at criticality, r — 0. Applying (1) and 
(2) to the scaling form (17) of the entropy, we obtain 







so 



V m p, 
spd 

1> 



$'(0)[ln(T o /T)]- d / ,/ ' +1 /(^) , (18) 
*(0)pn(T /T)]- , W- 1 . (19) 



Forming the ratio 0/c p , we find that the critical part of 
the Griineisen parameter behaves as 



(20) 



Equation (20) holds in the entire quantum critical region 
T > ToCxp[—\y x /r\ u ^] where the constant y x marks the 
crossover of $(j/). The sign of T does not follow from 
the scaling analysis, but because the entropy accumu- 
lates close to the finite-temperature phase boundary, we 
generally expect T > in the quantum critical region. 6 

Let us now analyze the scaling form (17) of the entropy 
in the WD quantum Griffiths phase, i.e., for r > and 
T < To exp[— (yx/rY^]- Using the large-argument limit 
of the scaling function Q(y), the density of surviving clus- 
ters is given by n(r, T) = Ar dv cxp[-cr^ ln(T /T)]. The 
resulting functional form of the entropy, 



S(r,T) = Ng(r)(T/T ) x ^ , 



(21) 



is identical to that found in (7) using heuristic rare region 
arguments, but the scaling analysis also gives A(r) = cr v ^ 
and g(r) = AsQT dl/ in terms of the distance to criticality 
and the critical exponents. Inserting g(r) and A(r) into 
(8), (9), and (10) leads to 







r = 



Asp 

V m pc 

As cr du+ ^{T/Tp) x ^ 
1 vi\> 
V m p-p 



r dv+v ^- l cvi){T/T Q ) x ^ \n(T /T) , (22) 

(23) 

ln(T /T) . (24) 



The prefactor of the logarithmic temperature dependence 
of r thus diverges as l/(p — p c ) at the QCP. 

As discussed above, the behavior on the ordered side of 
the transition, i.e., in the WO quantum Griffiths phase, 
is dimensionality and disorder dependent. Once these are 
fixed, the analysis can be performed in complete analogy 
to the WD quantum Griffiths phase. 

Conclusions — We have determined the Griineisen pa- 
rameter r at pressure-tuned QPTs in the presence of 
quenched disorder. At an infinite-randomness QCP, the 
critical contribution to V diverges as [ln(T /T)} 1+1 /^ 
with T — > 0. In the associated quantum Griffiths phases 
on both sides of the QCP, the rare region contribution 
to r behaves as ln(To/T) with a prefactor that diverges 
and changes sign at criticality (r < for p < p c and 
r > for p > p c ). Our results must be contrasted with 
the behavior at clean QCPs, where T diverges as a power 
of T at criticality but remains finite for all p ^ p c - 



4 



In many systems, a QPT can be induced by doping 
instead of pressure. If the main effect of doping is an ex- 
pansion or compression of the lattice, it acts as 'chemical 
pressure'. Close to criticality, the effects of pressure p 
and doping x can then be related via (p — p c ) = c(x — x c ) 
with c a constant. Defining the distance from criticality 
for such a transition as r = (x — x c )/x Cl this leads to the 
relation d/dp = (cx c )~ 1 d/dr. All our results for (3 and 
r thus hold if cx c is substituted for p c . 

If the transition is tuned by magnetic field instead 
of pressure, our calculations carry over to the magne- 
tocaloric effect Th defined in (4). In fact, by replacing 
VmPc by H c in the results for T, one obtains the cor- 
responding expressions for Th- Note, however, that our 
analysis assumes random- T c type disorder which does not 
break the order parameter symmetry. In magnetic-field 
tuned transitions in the presence of disorder, stronger 
random-field type effects may be generated. 30,31,32,33 
They would require a separate analysis. 

In our LGW theory, quantum Griffiths phases and 
infinite-randomness QCPs occur either for Ising symme- 
try without dissipation or for continuous 0(n) symmetry 
(n > 1) and Ohmic dissipation. 21 So far, we have fo- 
cused on the Ising case. The main difference for n > I 
is that ordered clusters act as (damped) quantum rotors 
rather than two-level systems. 29 In the Griffiths phase, 
this changes the prefactors g(r) and g(r) while the ex- 
ponents A(r) and A(r) remain the same. At criticality, 
the entropy (17) picks up an extra factor ln(Xb/T) from 



the entropy of a single rotor. It drops out in the ratio 
(3/cp. Thus, our results remain valid in the 0(n) case, at 
least at criticality and in the WD Griffiths phase. In the 
WO Griffiths phase, the rare region contributions to (3 
and c p will be overcome by conventional soft mode terms 
(because the dissipative 0(n) system is gapless). 

We now turn to experiment. CePdi-zRhi; features a 
ferromagnetic QPT that appears to be dominated by rare 
regions. 22,23 The phase boundary develops a tail charac- 
teristic of a smeared QPT, and at temperatures above 
the tail, magnetization, susceptibility and specific heat 
display nonuniversal power-laws as expected in a quan- 
tum Griffiths phase. Recent measurements of the ther- 
mal expansion 23 resulted in a very weakly temperature- 
dependent Gruneisen parameter close to the putative 
transition at x c « 0.87, in qualitative agreement with 
our theory. However, the variation of T with doping x 
differs considerably from our results. This may be caused 
by the fact that the doping is not isoelectronic. It thus 
not only acts as chemical pressure by inducing a lattice 
compression (as assumed in our discussion), but it also 
changes the electronic structure directly. To disentangle 
these effects one could prepare a sample with doping close 
to x c and then drive it through the QPT by pressure. 

We acknowledge discussion with M. Brando, M. 
Deppe, P. Gegenwart, and T. Westerkamp. This work 
has been supported in part by the NSF under grant 
no. DMR-0339147, by Research Corporation, and by the 
University of Missouri Research Board. 
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